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Abstract 

The impact of damping effect and external forces to the DNA breathing is 
investigated within the Peyrard-Bishop model. In in the continuum limit, 
the dynamics of the breathing of DNA is described by the forced-damped 
nonlinear Schrodinger equation and studied by means of variational method. 
The analytical solutions are obtained for special cases. It is shown that 
the breather propagation is decelerated in the presence of damping factor 
without the external force, while the envelope velocity and the amplitude 
increase significantly with the presence of external force. It is particularly 
found that the higher harmonic terms are enhanced when the periodic force 
is applied. It is finally argued that the external force accelerates the DNA 
breathing. 
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1. Introduction 



It is well known that the biological molecule functions such as transcrip- 
tion and replication can not be explained only by their static structure but 
also by their dynamical behavior [l|. The transcription processes begun with 
separation of the double helix into a single helix called denaturation process 
jij. The famous model for explaining thermal denaturation is the Peyrard- 
Bishop (PB) model [Sj] . The model consists of two chains connected by Morse 
potential representing the hydrogen (H) bonds. It has been shown that the H 
bonds stretching depends on the coupling constant and temperature, while 
in the continuum limit satisfies the Nonlinear Schrodinger equation. The 
solution of Nonlinear Schrodinger equation is implies for a self-focusing case 
could initiate the denaturation jsf. The intensive studies on PB model have 
been done for example, the molecular dynamics calculation j^, including 
anharmonics nearest-neighbor stacking interaction jsj, connection with con- 
formation of local denaturation jej, the effects of stacking interactions 0| and 
the Monte-Carlo simulation jsj. 

Even under physiological conditions, the DNA double-helix spontaneously 
denatures locally, opening up and fluctuating, and have a large amplitude 



localized excitations called DNA breathing [9|, llO|. It is well known that the 



thermal denaturation of double strands DNA depend on the solution that 



surrounding the DNA molecules Using this fact, it is necessary to take 
into consideration that the solving water does act as a viscous medium that 
could damp out DNA breathing. The impact of viscosity was investigated by 



Zdrakovic et al [12[. The behavior of DNA dynamics in viscous solution is 
described by the damped nonlinear Schrodinger equation. In the other hand, 
the experiment showed that the double stranded DNA can be separated 



by applying a fixed external force [13|. It was shown that the observed 



phase diagram for the unzipping of double stranded DNA is much richer 
than the earlier suggestion theoretical work. Therefore, it is necessary to 
improve the PB model to take into account a viscosity and the external forced 
simultaneously. To our knowledge, this problem has never been reported 
elsewhere. 

This paper discusses the viscous dissipation effect as well as an exter- 
nal force acting to DNA in the Peyrard-Bishop models. In the section ([2]) 
formulation of the impact of a damping and external force is derived. The 
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continuum approximation to describe DNA breathing is described in section 
([2D- Soliton solution based on the variational method is given in section (jl]). 
The breathing dynamics in term of solitary waves is discussed in section ([5]). 
The paper is ended with a summary. 



2. Peyrard- Bishop Model with damping and external forces 

Following the PB model, the motion of DNA molecules is represented by 
two degree of freedom Un and Vn which correspond to the displacement of the 
base pair from their equilibrium position along the direction of the hydrogen 
bonds connecting the two base in pair in two different strands j^. Making a 
transformation to the center of mass coordinate representing the in phase and 
out of phase transverse motions, X„ = (m„ + Vn) / V2 and ¥„, = [un — f„) / \/2 
respectively, the Hamiltonian of the PB model is given by , 

n n 

+^iYn^i-Ynr + §ie~^'^-~l)\ (1) 

where D and a are the depth and inverse width of the potential respectively. 
The momentum = MX^, Pn = MYn and k, is the spring constant. 

As mentioned above, the studies of PB models that included the viscosity 
was done by adding the term —e'jYn in the equation of motion (EOM) [l2| . In 
ref. |12| the nonlinear Schrodinger equation with viscous effect was solved to 
study the dynamics of DNA breathing. The interaction between the system 
with it's environment lead dissipation of energy. This means that the system 
is not longer conservative and reversible. The corresponding Hamiltonian 
formulation for dissipative system is called Caldirola-Kanai Hamiltonian in 



the form of a time-dependent Hamiltonian which defined as follow |14| . 
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H = e"^*^ + e^'V(x) , (2) 
2m 

where 7 = rj/M, rj is damping coefficient. The model has been used to study 
the quantum dissipation such as the quantization of an electromagnetic field 



inside a resonator filled by dielectric medium [15[, study of susceptibility 



for identical atoms subjected to an external force [16|, coherent states for 



the damped harmonic oscillator [16|, dissipative tunneling of the inverted 
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Caldirola-Kanai oscillator 17| and functional integral for non-Lagrangian 
systems 



In principle one can extend their approach to describe the denaturation 
processes in a dissipation system. We propose that the PB model with the 
damping effect and an external driving force F[t) is defined as follow, 

= J]e-^*^ + e^*^(X„+i-XO' + e^*F„(t)X„, (3) 

n 



n 

+ ^(e-f^"-l)2 + e^*F„(t)F., (4) 



Here -F„(t) is a conservative force. Substituting into the Hamilton equation 
yield, 

dPn M 



+ e^*^(e-t^"-l) + e-t^"+e^*F„(t). (5) 

Next, by substituting first equation into the second equation we find the 
corresponding EOM as follow. 



MYn + M^Yn = K{Yn+l-2Yn + Yr 



n—1 ) 



+ — (e-^^"-l)e-2>'"+F„(t). (6) 

In the mean time, the EOM for X„ satisfies the following damped harmonic 
oscillator with external forcing, 

MX„ + M7X„ - ft:(X„+i - 2X„ + X„_i) + = . (7) 

The damping term vn'^Yn is similar with pJi], where they add a new force 
= -^Yn in the EOM. 
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3. Continuum Limit Approximation 



The dynamical behavior of DNA breathing can be studied by applying 
the continuum approximation on the equation ([7]). We assume that the am- 
plitude of oscillation is small and the nucleotide oscillate around the bottom 
of the Morse potential but larg e enough due to nonlinear effect. We use the 
following approximation [H, [HI, 19 1, 



Yn = e^„: (e < 11 



(8) 



Substituting Eq.® into Eq.(|6]) and retained up to the third order of the 
Morse potential, we get 



(9) 



where Uq 



Cl 



2M ' 



3 

the average value of D and Fn = -F„(e\l/„). For a relatively long DNA chain, 
this equation can be simplified by taking full continuum limit approximation 
which should be valid as long as the solution under consideration changes 
rather slowly and smoothly along with DNA This approximation yields, 



IS 



^^^J 



dt 



+ C;t.(^ + + e ) + Fix, t) 



(10) 



where Cq = u^P and / is a length scale. The term of eF{x,t) and e'jd'^/dt 
can be treated as a perturbation by assuming that it contributes small enough 
to the whole DNA motion which should be dominated by the first and second 
terms yielding a solution of \1/ ~ ^tikx-uit)^ Then, using the multiple scale ex- 
pansion method, namely by expanding the associated equation into different 
scale and time spaces [20|, 



d_ 
dt 
d_ 

dx 



d d 



+ 



_d_ 

dxr 



+ e 



d 

dxi 



+ ■ 



(11) 
(12) 

(13) 
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and by substituting this expansion into Eq. (fTU]) . one obtains, 







dtl 







+ 2- 



-2Cl 



5^(1) 



r)2vI/(0) o 



^2^(0) ^2^(1) 



,^2^(0) 



ato^ti ° ax? 



+ 



(14) 



Since the first and second terms in the leading order (LO) and next to leading 
order (NLO) terms of Eq. fll4l) provide the harmonic solutions, while the 
remaining terms lead to non-harmonic solutions, then it is reasonable to 
consider, 



F = /(xi)e^('^^°-'^°*«) 



(15) 
(16) 
(17) 



Note that the time dependence on the right hand side of Eq. (ITTj) is intro- 
duced temporarily just for the sake of convenience, and there are additional 
charge conjugate terms in each equation. These lead to the order by order 
EOM, 



^(0) ^ 3 1^(1) |2 ^ 
i— — + Ai +A2- 



dr 



dr 



1^(1)1 



(18) 
(19) 

(20) 



where Ai = Co5/(2C^), A2 = 7/(2C^), A3 = 2D/C^,r = U = eto , 
^ = xi — {Cok/Cm)to, xi = exo and a dispersion relation Uq = D + Clk"^. The 
EOM of \E'''^-' is nothing else than the forced-damped nonlinear Schrodinger 
(FDNLS) equation. 
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4. Variational Methods 



In this section we discuss the solution of the FDNLS by means of varia- 
tional methods based on the Lagrangian formulation. It is well know that for 
the case with 7 = and / = 0, the Nonlinear Schrodinger equation admits 
the following traveling wave solution j^l, H, 19|. 



^(i)(^,r) = Aosech 



UeT) 



'i{ki;—UJT) 



(21) 



where, 



An 



(2A1A3) 



ut 



2Ai 
2Ai 



(22) 

(23) 

(24) 
(25) 



Here, Ug is the envelope wave velocity and Uc is the carrier wave velocity, 
satisfying ul — 2ueUc > 0. By using Eq.flTT]). Eq.(!T5l). Eq. (fT6l) we obtain the 
soliton solution as follow, 

'i/{x,t) =2¥^^cos{kx-uJot) + e\¥'^^\'^[3 + cos{2{kx-coot))] , (26) 

where is described by (pT]) . 

Based on the corresponding variational methods to solve the damped- 
forced nonhnear Schroedinger equation, one can use the solution (12T!) as the 
related basic form and considering its amplitude, width, phase velocity and 
the position of the soliton to be time dependent 20|, l2ll ? ]. For convenient, 
let us write the 1-soliton in the following form, 

vl/«(e,r) = r/(r)sech[r7(e + C(r))]exp(-^[^(r)e + 0(r)]) . (27) 

The dynamics of rj, 6, ( and (p function can be obtained by using the vari- 
ational methods. The Lagrangian that satisfy the Euler-Lagrange equation 



20 




dl 



c}v[/(i)* 



. 



(28) 
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for the FDNLS is given by, 

I = _ _ I |2 +A3 I |4 

+ _ ^a)*,j,(i)) _ _^ (29) 

Substituting equation fl27j) into the equation fl29|) and using L = /ci^, 
sech(a^)(i^ = vr/a and sech^(a,^) tanh(a^)(i,^ = yield, 

4 

L = 2r]Ce + 7]^ + 2iA2VCd + iA2V<j) + 4Air/^2 + -A^r]^ - 7]F , (30) 

3 

where 

F= r ('/V[^W«+*(^)1 + f*e~'W-)^+'f'(-A sech[r/(e - C{rM. (31) 

J —CO 

Eq. ( I30i) is the Lagrange function in term of 9, rj, and (. The EOM can 
be easily obtained by solving the Euler-Lagrange equation, 

d f dL\ dL 



MUX) OX »• 



where X = (77, 6, (j), (). Substituting the Lagrange function given by Eq. fl30l) 
into Eq. fl32l) leads to the following, 

dF 

2(l + zA2)r7C + 2(l + ^A2)C?7 = 8Airi9 - r]— (33) 

8F 

2(1 + ^A2)C^ + (1 + «A2)0 = A{A,e^-Asv')-V-Q^-F (34) 

OF 

{l + zA2)r] = -r^— (35) 

2(1 + 2A2)^ = (36) 

Further, by substituting Eq.( l35|l and Eq. (l36|l into Eq.( l33|l and Eq. (lM|l re- 
spectively, we find 

dF dF 

{l + iA,)^C~^A,e = 2r/— -— (37) 
(1 + ^A2)0 - 4(Aie2 - Aar^^) = - - ^ • (38) 
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Generally, the analytic solution of the equation can not be obtained, 

but it is still possible to find it for special condition. In the next section we 
restrict ourselves to special case of /. 

First, let us consider very special case in which the damping factor and 
external force are vanish (A2 = and F = 0). As a result , Eg. (1351) show 
that the function of 17 = r^o is a constant and Eq. flSB]) yield 9 = 9q which is 
also a constant. Eq. flS^ and Eq. fl51|) lead to a simple solutions, 

C = 4Ai^or (39) 
= A{K,el - K^riDr . (40) 

Such that single soliton have the form of, 

r) = r/osech [r/o(e + 4Ai^or)] exp (-^[^0^ + K^idl - ^^vDA) • 

(41) 

Clearly, this is a single soliton solution of the conventional Nonlinear Schrodinger 
equation. 

In the second case, let us ignore the external forces (F = 0). Again, as 
a result, Eq. fl35|) show that the function of r/ = r^o is a constant and Eq. fl36l) 
yield 9 = 9o, constant as well. The time dependent variables are ( that 
represent of the velocity of soliton and (f) the phase of soliton. The equation 
of motion have a simple form as follow, 

{I + tA^)^ - AAi9o = (42) 



This yield. 



:i + ,A,)2 - 4 {A,0l - A„i} = (43) 



,M . 4(^^,1 A,), (45) 



*-2 

Finally, the single soliton solution is. 



X exp (-(^Aar) exp (-^(^0^ + 0r)) , (46) 
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where C = 4Ai^o/(l + A^) and = A{K^el - K^r]l)/{1 + A^). 

For the third case, we consider another simple case namely when F = Fq 
is a constant. Again we get r] = 7]q and 9 = 9q are constants. The EOMs in 
Eq.f l37|) and Eq.f pSj) have simple forms as follow, 



This yields. 



CO 



(1 + 2A2)^-4AA 
dr 



2 







4 (Aigg - Aar/g) 
(1 + A2^ 



1 + A2^ 



^2; V-L-r^i-2 

Finally, the single soliton solution for this case is, 
vl/«(e,r) 



(47) 
(48) 

(49) 
(50) 



X 



r/oexp (-0A2r) exp (F*r) 
sech[r7o(e + C^) + tr]oC^2r] exp [-^(^0^ + 0^ - (51) 



where F* = Fq/{1 + Ag). The external force express in the positive 
exponential term, this show that it will increase the amplitude of the soliton. 

In the fourth case, we assume that the amplitude and the phase is a 
constant ij] = rjo) and {6 = 9q), then the EOM becomes. 



(1 + iA2)voC - 4Ai^o 







dF 



(l + ^A2)0-4(Ai0o-A3%') = C^-i^ 



(52) 
(53) 



Further, let us assume that the external force is expressed specifically by 

f = foexp{-z9oO , (54) 
with /o is a constant. Substituting into Eg. (1311) yields. 



{foe 



-ieo5„i[0o€+0(r)] 



+ /oe 



ieo(^-im+HT)] 



) sech[r/(e - CirM 



(A\ 

COS(0j 

71 



(55) 
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Such that we find, 



1 - V« „osW + l(Mz/^. (57) 



7r(l + iA2) ' ' (l + iA2) 

The Eq. fl57|) will be solved numerically. It is important to note that the 
solution of Eq. fl57|) is a complex function i.e. (p = (pR + such that the 
soliton profile is given by, 

^^'He, r) = vo exp {(Pit) sech[r]o{^ - (r) + tVoC^2r] exp [-1(60^ + (Prt)] , 

(58) 

5. Nonlinear Dynamics of DNA Breathing 

The external effects of DNA usually give inhomogeneities in the DNA 
model. The inhomogeneity in stacking energy is found to modulate the width 
and speed of the soliton which depend on the nature of the inhomogeneity 
[? ]. The author used the dynamic plane-base rotator model by considering 
angular rotation of bases in a plane normal to the helical axis. They found 
that the DNA dynamics is governed by a perturbed sine-Gordon equation. 
In this paper we found that the inhomogeneities of the DNA breathing dy- 
namics governed by the forced-damped Nonlinear Schrodinger equation. The 
solutions of the homogeneous case represent a large amplitude with localized 
oscillatory mode appears as a good explanation of the breathing of DNA and 
must be spontaneously formed [H, 0, 10 . 



Let us discuss the previously considered four simple cases. In the first 
case, it's shown by using transformation tjq = Aq and 6q = l/2ue/Ai that 
the solution is just the Eq.( !2T]) . We simulate the solution for the model 
parameter given by k = 8Nm, M = 5.1 x 10~ ^^fcq , a = 2 x 10^°m~^, 
D = O.leV and / = 3.4 x 10~^°m is length scale 12]. The system of unit 



{A°,eV) defines a time unit (t.u.) equal to 1.021 x 10~^^s [22|. The solution 
demonstrate a sort of a modulated solitonic wave where the hyperbolic and 
cosine terms correspond to the wave number of the envelope and the carrier 
wave respectively. 

In second case with F = and damping constant 7 = 0.05kg/ s, the be- 
havior of the breathing is depicted in Fig.([T]). The figure show that the the 
breathing propagation along the DNA molecule is effected by the damping. 
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Figure 1: The DNA breathing in the second case where x denotes the continuum base pair 
in the present model (black) and the original PB model (red) with g — 0.05 and F = 0. 



It is shown from the figure that the damping term decelerates the propa- 
gating sohton while retaining its amplitude profile. This indicates that the 
corresponding damping term does not affects the soliton mass. 

Now, consider the solution with the present of the damping effect and 
F ^ where the solution is depicted in Fig.(j2]). The figure is generated 
with F = 15pN and a damping factor 7 = 0.05kg/ s. The external force 
tend to increase the breathing amplitude and damped out by the damping 
effect. The consequences of this matter, the present of external force tend to 
increase the envelope velocity significantly. The velocity is increase around 
Aue = AC, I Ar ~ lorder in the present of external force F = 15pN. 
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t=0.02 




Figure 2: The solitonic solution of DNA breathing in the model with damping effect and 
external force F (black) and in the original PB model (red) with 7 = 0.05 and F = 15. 



Our result is supported by the report given in ref. 
cell, DNA strands are separated by the external force 23 




results. In a 
or in chemical 

terms, by enzymes whose interactions with DNA make strands separation 
thermodynamically favorable at ambient temperature ^5|. They showed that 
the two strands of double-stranded DNA can be separated (unzipped) by the 
application of 15pN force applied at room temperature. Their model predicts 
that the melting temperature should be a decreasing function of applied force. 
The paper show that the external force can increase the amplitude of the 
breathing and may separate the double helix into single helix. 

Finally, let us consider the fourth case. The behavior of this case described 
by on the numerical solution of Eg. fl57p which gives a complex function of 
(j). The real part of the solution is associated with the carrier wave, while 
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0.2 



0.4 



0.6 



0.8 



T/t„ 



Figure 3: The solution of Eq. (j57p where the real and imaginary parts are shown by the 
dot and solid lines. 



the imaginary part with the envelope velocity respectively. The profile of the 
real part and the imaginary part can be seen in Fig.(|2]). The real part profile 
seems to have an anti-sigmoid-like function. In the mean time, the profile 
of the imaginary part have negative value, which means that the solution 
propagates to the left direction. 

At this point we can consider the solution of FDNLS equation by using 
Eq.f lFTI) and this is depicted in Fig.(|l]). Soliton solution of FDNLS show an 
increasing of the amplitude and its velocity generally tend to increase with 
time. The amplitude is relatively high and the shape of soliton tend to be 
step form during it's propagation. 

Now, the DNA breathing of the fourth case using the solution of "^^^^ given 
in Fig. dl]) is depicted in Fig.^S]). The amplitude of the DNA breathing is 



around 2nm where the result is similar with ref. [12| for the same parameters. 
The amplitude of the DNA breathing is relatively increase when its propagate 
to the left direction. The high amplitude show that the DNA tend to be 
unzipped and finally completely separate into single helix. 

Finally, it is interesting to note the case if we decrease the envelope and 
carrier velocity respectively in one order. The results is depicted in Fig.d^]). 
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Figure 4: The solution of soliton in the FDNLS (black) and the NLS (red) with 7 = 0.05, 
F = 15e*^, Ve = 10^ ni/s and Uc = 4 x 10"* ni/s. 



In this case, the solution of FDNLS propagates more slower than the previ- 
ous one. At r = 0.3 the soliton with periodic external force tend to increase 
its amplitude and velocity. Further at r = 0.6 the amplitude increase signif- 
icantly and the corresponding form is change significantly. At r = 0.7 the 
higher harmonic term is developed and completely forms at r = 0.9. As the 
times goes up the amplitude tend to decrease and disperse into wider form 
and then the higher harmonic term is generated and increase the amplitude. 
It is interesting to pointed out that the harmonic term come from the the 
solution of Eq. flFTl) that is a nonlinear equation. 

The DNA breathing corresponds to the "^^^^ above is depicted in Fig. ([7]). 
The result shows that the amplitude is about 0.2nm which is the same with 
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Figure 5: The DNA breathing propagation for the fourth case in the present model (black) 
and the original PB model (red) with 7 = 0.05, F = 15e'^, Ve = 10^ m/s and Mc = 4 x 10^ 
m/s. 



the result found in ref. 22|] in which they used the Forinash-Cretary-Peyrard 
model with helicosity taken into account. They showed that the opening of 
the double helix of the DNA itself is controlled by the resonance mode. 

As the result shows that the periodic external force and the damping effect 
generate a higher harmonic term in the dynamics of FDNLS solution then 
it can be concluded that this phenomenon is responsible for the dynamics 
of the breathing of the DNA. At this point, we can see that in the early 
propagation process the DNA breathing tend to decrease its amphtude and 
disperse into a wider form. The condition changes when the higher harmonic 
term of FDNLS soliton begin to develop which leads the amplitude of the 
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Figure 6: The solution of FDNLS (black) and NLS (red) with 7 = 0.05, F = 15e'«, K ^ 4 
Nm, Ve ~ 10'* m/s and Uc = 4 x lO'^ m/s. 



DNA breathing increases significantly. 
6. Summary 

The impact of viscous fluid and external forces to the Peyrard-Bishop 
DNA breathing is investigated. We have proposed a PB model with the 
damping effect and the external driving force which is described by an ex- 
tended time-dependent Caldirola-Kanai Hamiltonian. Taking full continuum 
approximation and using the multiple scale expansion method, the EOM is 
nothing else than the FDNLS equation. Assuming small perturbation of 
damping and external forces, the FDNLS equation can be solved using the 
variational methods. The analytical solution have been obtained for only spe- 
cial cases. In the case with damping factor and without external force, the 
breathing propagation is decelerated by the damping effect. In the presence 
of external force, the velocity and the amplitude increase significantly. It is 
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also found that the higher harmonic terms are enhanced when the periodic 
force is apphed. 

These results shows that the external force contributes constructively to 
accelerate the separation of double helix into single helix. More comprehen- 
sive numerical investigation of the variational equation found in this paper 
is still in progress. 
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